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Probing the neutron star spin evolution in the young SMC Be/X-ray 
binary SXP 1062 
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ABSTRACT 

The newly discovered Be/X-ray binary in the Small Magellanic Cloud, SXP 1062, provides 
the first example of a robust association with a supernova remnant (SNR). The short age 
estimated for the SNR qualifies SXP 1062 as the youngest known source in its class, r ~ 
10 4 yr. As such, it allows to test current models of magneto-rotational evolution of neutron 
stars in a still unexplored regime. Here we discuss possible evolutionary scenarios for SXP 
1062 in the attempt to reconcile its long spin period, P = 1062 s, and short age. Although 
several options can be considered, like an anomalously long initial period or the presence of a 
fossil disc, our results indicate that SXP 1062 may host a neutron star born with a large initial 
magnetic field, typically in excess of ~ 10 14 G, which then decayed to ~ 10 13 G. 
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1 INTRODUCTION 

Be/X-ray binaries (or BeXBs for short) form a subclass of the high- 
mass X-ray binaries (HMXBs) in which the neutron star (NS) com- 
panion is a Be star, a spectral class B giant/subgiant with emission 
lines and large IR flux. These peculiar properties are explained in 
terms of an equatorial disc, formed by matter lost by the rapidly 
rotating Be star. X-ray emission is believed to be powered b y ac- 
cretio n of material in the equatorial disc onto the NS (see e.g. lReid 
1201 lL for a recent review). 

BeXBs are both transient and persistent X-ray sources. Tran- 
sient systems are characterized by type I-II outbursts during which 
their flux increases by a factor 10-10 4 over the quiescent level. 
They typically contain a not-too-slow NS (P < 100 s) on a mod- 
erately eccentric orbit, P OT h i$ 100 d, e > 0.2. On the other hand, 
persistent BeXBs exhibit a rather flat lightcurve, lower X-ray lu- 
minosity (L ss 10 34 -10 35 erg s _1 ), l onger spin a nd orbital periods 
(P > 200 s, P orb > 200 d; see again lReisfeOllI) . 

There are presently about 30 well-established BeXBs in the 
Galaxy, plus ~ 40 candidates. In addition, ~ 50 sources (plus ~ 20 
candidates) ar e known in the Smal l Magellanic Cloud (SMcfl 
Very recently iHenault-Brunet et al.1 d2012l) reported the discov- 
ery of a new BeXB in the Wing of the SMC. The new source 
(SXP 1062) has the typical properties of a persistent BeXB: a B0- 



* E-mail: sergepo lar@ gmail.co m 

1 Figures from iReid 1201 ll) and the BeXB o nline catalogue at 
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0.5(III)e+ donor, L ~ 6 x 10 35 erg s" 1 , P ~ 1062 s, and an 
orbital period P or b ~ 30 d, as estimated from the Corbet di- 
agram icorbetetai1l2009l) . What makes SXP 1062 to stand out 
amongst its kinship is its (likely) association with a supernova rem- 
nant (SNR). BeXB-SNR associations have been already reported 
(again in the SMC) but i n all previous cases they a ppear uncertain 
dHughes & Smithll 19941 ; Iboe. Haigh & Reidl2000l) . In the case of 
SXP 1062 the association looks robust and allows for the first time 
to estimate the NS a ge in a BeXB from that of the parent SNR, 
r ~ 2-4 x 10 4 vr iHenault-Brunet et alj|2012l) . The suggested 
association of SXP 1062 with a SNR has been further strength- 
ened by a reanalysis of the same XMM-New ton datasets, supple - 
mented with optical and radio observations, by |Haberletai] ( l2012l) . 
Their estimate for the SNR age, 1.6 x 10 4 yr, is even shorter 
than, a lbeit fully compatib le with, that of IHenault-Brunet et al.l 
d2012l) . Eaberl et all J20I2I) were also able to measure the source 
period evolution, which results in a positive (i.e. spin-down) rate 
P~3x lO^ss- 1 ~ 95syr-\ 

The long spin periods (P > 1000 s) of some persistent BeXBs 
have been for a long time a major issue. According to the stan- 
dard picture, there are four stages in the spin evolution of a neu- 
tron star embed ded in a medium : ejector, propeller, accretor and 
georotato^f] (e.g. lLipunovlll992l) . Once the NS entered the accre- 
tor stage after a short propeller phase, its spin period quickly set- 
tles at an equilibrium value, P eq . In the conventional model, based 
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on lDavies & Pringlel d 198 lh results, the star dipole field should be 
B > 10 14 G to have P eq > 10 3 s, unless the accretion rate is or- 
ders of magnitude below what needed to account for the observed 
X-ray luminosity. Since observations support the presence of a NS 
with standard magnetic field (at least in some BeXBs), the sub- 
sonic propeller stage, which ca n delay the onset of accretion until a 
much longer period is reached (ikhsanov 2007j), has been invoked 
to explain ultra-slow BeXBs (see Reidl2011 ). More recent investi- 
gations of wind-fed accretion onto magnetized NSs indicate, how- 
ever, that the equilibrium period can be as high as ~ 1000 s even 
for 5 m 10 12 -10 13 G and M « 10 16 g s -1 , as expected in BeXBs 
JShakura et alj2012h . 

Whatever the details of the braking torques, the previous ar- 
gument implicitly relies on the assumption that the present age 
of the source is long enough for the NS to have entered the pro- 
peller stage. Unless the accretion rate is way above typical, a NS 
in a BeXB with B > 10 12 G starts its evolution in the ejec- 
tor (or pulsar) phase. Its duration can be roughly estimated as 
Tej > 10 6 (B/10 12 G)" 1 (Af/10 15 gs- 1 )" 1/2 yr(seeSection[2l(. 
This is comfortably below (by a factor ~ 10) the lifetime of the 
Be companion, so there is ample room for the binary to start an 
accretion-powered X-ray stage. In the case of SXP 1062, however, 
it would be impossible for the NS to enter the propeller stage (and 
hence to become an accretor) in a time as short as a few x 10 4 yr, 
the estimated SNR age, for typical values of B and M. The accre- 
tion rate in SXP 1062 is M = L/rjc 2 ~ 6 x 10 ls gs _1 for an 
efficiency r\ = 0.1, so this points to a highly magnetized NS, with 
an initial magnetic field substantially above 10 12 G. 



2 SPIN EVOLUTION IN SXP 1062 
2.1 Spin-down torques 

In th e standard picture of NS spin evolution (see e.g. iLipunovl 
1 19921) . the transitions among the different stages are regulated 
by the relative values of some characteristic radii. In the ejector 
phase the light cylinder radius, Ri = cP/2-k ~ 5 x 10 9 P cm, 
is typically smaller than the gravitational capture radius, Rg = 
2GM/V 2 ~4x 10 11 V^oo cm, where V is the velocity (and p the 
density) of matter far from the staff A NS with M N s = 1.4 M , 
Rns = 10 km and moment of inertia / = 10 45 gem 2 is assumed 
henceforth. The transition to the propeller stage occurs when the 
ram pressure, Pd yn = pV 2 /2, balances the outgoing flux of elec- 
tromagnetic waves and relativistic particles, Ppsr = E /(A-kR 2 c), 
at Rg (E is the rotational energy loss rate of the pulsar). Once mat- 
ter crosses Rg it is gravitationally captured and quickly reaches the 
light cylinder radius switching off pulsar emission, since the dy- 
namical pressure (assuming a nearly spherical flow) rises as R~ 5 ^ 2 , 
while the relativistic momentum flux goes like R~ . 

The critical period for the transition follows by requiring that 
Pdyn(Rc) = Ppsr(Rg) together with the standard expression 
for magneto-dipole losses, E = 8tv 4 B 2 R% s sin 2 a/(3c 3 P 4 ) and 
mass conservation, 4tvR g pV = M, 
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3 Hereafter V 3 00 = (V/300 km s" 1 ), P orb300 = (P ort /300 d), 
P1000 = (P/lOOOs) and Nx is used for the quantity N in units of 
10 x cgs. 



since the NS is not accreting, here M is more properly defined as 
the matter flow rate in the surrounding medium and we assumed 
a nearly orthogonal rotator (sin a ~ 1). From the magneto-dipole 
formula, assuming constant B and very short initial period, the time 
it takes the NS to spin down to P ej , i.e. the duration of the ejector 
stage, can be evaluated 



3/c 3 P 2 



ej l6n 2 B 2 R% s 



1.5 m: 



l/2 T/ -l/2 



B^i Myr. (2) 



The dipole fiel d in a wind-fed NS has been estimated by 
Shak ura et al. U20T1 under the as sumption that the sta r is spinning 
at the equilibrium period (see also lPostnov et alj|201ll) . 
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The previous expressions give r e j ~ 0.2 Myr for SXP 1062, about 
an order of magnitude larger than the estimated SNR age. Our ref- 
erence value, V ~ 300 km s , is a conservative estimate for the 
typical velocity of matter ejected from hot stars in binaries (see 
iRaguzova & Lipunovll 19981) . Actually, since Eq. l[3} depends rather 
strongly on V, higher velocities would result in a lower magnetic 
field and longer r e j . 

Within this framework, an obvious possibility to shorten the 
ejector phase in SXP 1062 is to invoke a higher dipole field. How- 
ever, if the present field is that given by Eq. $3$ this implies that 
B must have been stronger in the past and then decayed to its 
present value. The issue of magnetic field evolution in isolated NSs 
and its observable consequences has been recently addressed in de- 
tail, especially in connection with strongly magnetized objects, or 
magnetars (eg.lPons. Miralles & Geppertll2"009t|Popov et alj201ol: 



iTurolla et al.ll201ll and references therein). In these studies it is 
assumed that magnetic field decay takes place in the stellar crust 
and is driven by Hall/Ohmic diffusion. Both processes are strongly 
density- and temperature- dependent, so magnetic and thermal evo- 
lution are necessa ry coupled, and a multidimension al numerical ap- 
proach is needed jPons. Miralles & Gepperlj[2009h . 

Since we are interested in tracing the spin evolution of SXP 
1062 prior it entered the accretor stage, we can treat the NS as iso- 
lated for the sake of magnetic evolution. In order to avoid com- 
plex numerica l simu lations, we adopt the simplified approach of 
lAguilera et al .1 d2008l) in which the dipole magnetic field decay is 
described by the analytical law 



B(t) 



B exp (— t/ro) 



1 + (to/t h )[1 - exp (-t/ro)] 



+ B 
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where Bo is the initial field, Bfi n is the relic field, to and th are 
the characteristic timescales of Ohmic and Hall decay, respectively. 
The period evolution in the ejector stage is given by the (general- 
ized) magneto-dipole formula 



87r 2 £ 2 i4 s (l + sin 2 a) 



(5) 



3/c 3 P 

( Spitkovsky 2006), where B is given by Eq. ©■ 

The NS enters the propeller phase as soon as the dynamical 
pressure exerted by the incoming material overwhelms the pulsar 
momentum flux at the gravitational radius, as discussed above. If 
no stable equilibrium exists, matter will reach the light cylinder ra- 
dius on a free-fall timescale, and proceed inwards if the ram pres- 
sure overcomes the magnetic pressure of the dipole field, P mag = 
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B /8tt. lEksi & Alpail J2005h have shown that stable matter con- 
figurations can be present even if the inner boundary of the flow, 
Rin, is outside Ri. This occurs when Ri„ < R cr u = F[o)Ra, 
with F(a) ^ la function of the magnetic tilt angle anqj Ra = 
(1/2)B 4/7 R 12/ s 7 /(8GM NS M 2 ) 1/7 .For R t < R ln < R crlU the 
pulsar is still active and the st ar is spun down by m agneto-dipole 
torque. Under these conditions lEksi & Alpaj d2005l) provide an ex- 
pression for Ri n in terms of a and Ra- When i?, n < Ri, instead, 
the flow is truncated at the Alfven radius, where Pd yn = Pmag, the 
pulsar is turned off and only propeller torques act. 

The propeller physics is very complicated and no unanimous 
consensus exists about the form of the torque. However, since spin- 
down is expected to be very efficient in the propeller phase, its du- 
ration is quite short. The period rapidly increases and the corotation 
radius, R co = [GM NS / '(4tt 2 )] 1/3 P 2/:i , moves quickly outwards 
until it matches the Alfven radius, at which point accretion begins 
and the period freezes at the equilibrium value. While the latter de- 
pends on the propeller/accretion mechanisms, and several options 
have been discussed (see Sec. [I}, the duration of the ejector phase 
does not. For this reaso n we consider here only the "maximally 
efficient torque" (see e.g. iFrancischelli & Wiierskdo^ . as an illus- 
trative example: 

p2 

P = -MR 2 n [Q K {Ri n )-2Tv/P]— (6) 

TVl 

where £Ik is the keplerian angular velocity. 

Finally, to follow the period evolution in the accretor stage, 
Ri n < R co , we assume the se ttling ac cretion regime recen tly pro- 
posed bv lShakura et al.l d2012l , see also lPostnov et al .11201 ll) 

P = - f4Mf 2 " )/u - CMU 11 ] ^ ■ (7) 

L J 27Ti 

Here A ~ 2.2 x W 32 K^B^Rns e)^ 11 Vf ^P-t 30 o> C ~ 
5.4 x 10 31 A'i(Bi2-Rjvs6) 13/11 -Piooo ( DOth in C 8 S ™its); the con- 
stant A'i and the index n were fixed to 40 and 2, respectively. The 
accretion torque can produce either spin-up or spin-down and it is 
interesting to note that in the spin-down phase P increases expo- 
nentially, with a typical timescale ~ lOOB^ 13 '' 11 M^^ 11 yr. 

2.2 Results 

We solved numerically the equation for the period evolution in the 
three stages (Eqs. (5), Qj, |7J) starting from to = 0.01 yr with an 
initial period Po = 0.01 s. The accretion rate was fixed to the value 
derived from current observations, Al — 6 x 10 15 g/s, together 
withPorb = 300 d, V = 300 km s~\ to = 10 6 yr and B fin = 
8xl0 12 G. Several cases were then computed varying the magnetic 
dipole angle a, the Hall timescale th and especially the initial field 
Bo - Figure[T]illustrates the results for a = 10°, th = 10 3 yr and 
B = 4xl0 14 , 10 14 ,7xl0 13 ,4xl0 13 , 10 13 G. A common feature 
to all evolutionary tracks is a very rapid propeller stage which is 
followed by an even more rapid spin-down phase as the NS enters 
the accretion regime. The period then saturates at its equilibrium 
value. The decrease in P seen at later times for the larger fields is 
due to the dependence of P eq on B and to the fact that the field 

4 Although Ra formally coincides with the Alfven radius, we remark 
that the latter exists only inside the light cylinder radius. Nonetheless, the 
present definition is correct, provided that Ra is not associated to the equi- 
librium of dynamical and magnetic pressure. 




time (yr) 

Figure 1. The spin period evolution for a = 10°, th = 10 3 yr and 
S = 4 X 10 14 , 10 14 . 7 X 10 13 , 4 X 10 13 , 10 13 G (solid lines, from 
top to bottom). The shaded areas mark the age (vertical strip) and period 
(horizontal strip) of SXP 1062 with the respective uncertainties. 

is still decaying (see Sec. 12. It . The time variation of B is also 
responsible for the deviations of P(t) from a power-law in the later 
ejector phase, which are, again, more prominent for large B. 

The main information Figure [JJ conveys is that a quite large 
initial field is required in order for SXP 1062 to enter the pro- 
peller phase (and quickly start accreting) in a time as short as a few 
x 10 4 yr. For the case shown here it has to be Bo > 10 14 G for this 
to occur. This result is not very sensitive to the actual choice of a 
and th - Both increasing a and decreasing th results in a somehow 
higher value of the minimum initial field, which is however in all 
cases in the range ~ 10 14 -4 x 10 14 G. The conclusion that SXP 
1062 harbours an initially strongly magnetized NS seems therefore 
quite robust. 

We stress that our main goal is not to reproduce within the cur- 
rent model the observed value of P at the present age. Although, for 
completeness, we followed the spin evolution also in the propeller 
and accretor stages, the treatment we employed is necessary ap- 
proximated and contains some degrees of arbitrariness in the choice 
of the propeller/accretor torques. Our computed evolution past the 
ejector phase has mainly an illustrative purpose and has to be taken 
with caution. 



3 DISCUSSION 

The quite short age implied by the association of the newly 
discovered BeXB in the SMC SXP 106 2 with a SNR 
l lHenault-Brunet et al.l [20121; lHaberl et at] 1201 2|) raises a number 
of questions about the properties of the neutron star and its evo- 
lutionary status. Normally, one expects that accreting X-ray pul- 
sars spin close to their equilibrium period. H owever, for a such a 
young system this appears far from granted. lHaberl et al.l J2012h 
reported a large spin-down rate for SXP 1062, which may suggest 
that P eq has not b een reached yet . According to the standard evo- 
lutionary scenario dLipunov|[l992l) . the (maximum) spin-down rate 
in the accretor stage is P ~ 2nB 2 R% s /(GMI) which implies 
B w 3 x 10 14 G for P » 100 syr~\ On the other hand, if Eq. 
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(0 is used to estimate the magnetic field, a much lower value is 
obtained, B ~ 10 13 G, very close to what is predicted assuming 
that the source spins at the equilibrium period (see Eq. (3)). This 
supports a picture in which the NS actually rotates close to P eq . 
A further argument in favor of this is the very short duration of 
the spin-down phase in the accretor stage, P/P < 100 yr, which 
makes it very unlikely to catch the source in this state. Our con- 
clusion is that both the young age and the large spin-down rate of 
SXP 1062 argue in favor of an initially highly magnetic NS which 
experienced field decay. 

Details of evolution past the ejector stage are uncertain, so a 
fine tuning of the model discussed in Sec.[2]is entirely premature. 
Still, some general considerations can be made. The main one is 
that the timescale required for reaching P eq after the NS left the 
ejector stage is very short. A further point to be stressed is that the 
zero, or negative, period derivative at P = P eq (see FigureQJ does 
not account for variations in M due to orbital motion in the BeXB 
and irregularities in the wind. This can result in quite large values of 
P and rapid changes in the period derivati ve which are indeed m ea- 
sured in wind-fed X-ray binaries (see e.g. lBildstein et al.ll 19971) . In 
this respect the spin-down rate in SXP 1062 is large but not excep- 
tional ly so. For examp l e, GX 301-2 i s know n to have a larger value 
of P dKoh et aljri997h . lHaberl et af] j2012h measured the average 
period derivative over a time-span less than a month, ~ one tenth 
of the orbital period. Fluctuations in P can occur on a timescale of 
days, so a more complete dataset is definitely needed to assess the 
real nature of the period variation. 

Alternative scenarios to explain the long period and short age 
of SXP 1062 can be envisaged. For instance. lHaberl et al.l d2012h 
suggested that the NS could have been born with an initial period 
much longer than ~ 0.01 s. The value of Po can be evaluated by 
requiring that the transition period given by Eq. l[T) is reached in 
less than the source age, and it turns out to be ~ Is for the P-field 
corresponding to P eq [Eq. {3)]. If this is the case no field decay is 
required. Although possible, no compelling observational evidence 
for such long initial periods in NSs exists. There are hints that some 
central compact objects (CCOs) in SNRs may have a current pe- 
riod very close to the initial one, Po ~ 0.1 s. These sources, how- 
eve r, are suspected to host we akly magnetized NSs (B ~ 10 11 G; 
e.g. lGotthelf & Halpernll201Cl) . Whether the low field and the long 
period are related, depending on the conditions under which the 
NS formed, is still unclear. Still, the estimated value of B in SXP 
1062 is much higher and the initial period required in the present 
case is more than one order of mag nitude longer. Very recently 
iKnigge, Coe & Podsiadlowskil d2012l) presented evidence that the 
population of BeBXs is bimodal, with the two sub-populations hav- 
ing distinct typical spin and orbital periods, and eccentricities. They 
suggest that the NSs in the two groups are produced through differ- 
ent channels, iron-core-collapse and electron capture supernovae, 
with the long spin period sources associated to the former channel. 
Whether this may result in a population of long Po NSs is a still 
open question. 

Another possibility is that the NS in SXP 1062 could, at least 
in the early stages following its formation, have been surrounded 
by a debris (or fossil) disc, left by the parent supernova explosion. 
The existence of such a disc could also lead to rapid spin-down 
and large period. This, in addition to a super-strong field, was sug- 
gested to explain t he ultra-long period ( ~ 6.67 h r) in the enigmatic 
source RCW 103 fc)e Luca et al J 200^11120071) . Although this re- 
mains a possibility worth of further i nvestig ations, p reliminary cal- 
culati ons, based on the model byjEJ J2007I see also lEsposito et all 
1201 lh , indicate that a large initial field (B > 10 14 G) is still re- 



quired, unless the disc is quite massive. For Mdisc (0) ~ 10 2 Mq , 
the NS can enter the propeller stage earlier than ~ 10 4 yr also for 
B ~ 10 13 G. 

If SXP 1062 indeed contains an initially strongly magnetized 
neutron star, then studies of this system can shed light on the ori- 
gin o f magnetars. In the standard scenario jDuncan & Thompson! 
Il992h super-strong fields are generated via dynamo processes. This 
requires rapid rotation of the proto neutron star. Primary compo- 
nents in high-mass binary systems are n ot expected to form rapidly 
rotating cores at the end of their lives l IPopov & Prokhorovll2006l : 
iBogomazov & Popovll2009T) . Up to now no strongly magnetized 
comp act objects have been identi fied in binary systems with cer- 
tainty. IChashkina & Popovl J20l3) have recently derived estimates 
of the B-field in HMXBs using IShakura et all J2012h formula 
and no evidence of ultra-high fields was found. It has bee n sug 
geste d that supergiant fast X-ray transient s (SFXTs, see e.g. Sidoli 



2011 , for a review) m ay host a magnetar 
20081) . More recently, (Torres et all J201 



JlBozzo. Falanga & Stella 
2i) reported magnetar-like 



behaviour from the peculiar binary LS I +61 303. However, no defi- 
nite confirmation has been given yet, and, in this respect, SXP 1062 
may be a unique example. The existence of a (evolved) magnetar in 
a high-mass binary system will pose new challenges on the origin 
of such neutron stars. 
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